P-matrices in orbit spaces and invariant theory 



V Talamini 

INFN, Sezione di Padova, via Marzolo 8, 1-35131 Padova, Italy 
E-mail: vittor ino . talaminiOpd . inf n . it 

Abstract. In many physical problems or applications one has to study func- 
tions that are invariant under the action of a symmetry group G and this is best 
done in the orbit space of G if one knows the equations and inequalities defining 
the orbit space and its strata. It is reviewed how the P-matrix is defined in terms 
of an integrity basis and how it can be used to determine the equations and in- 
equalities defining the orbit space and its strata. It is shown that the P-matrix is 
a useful tool of constructive invariant theory, in fact, when the integrity basis is 
only partially known, calculating the P-matrix elements, one is able to determine 
the integrity basis completely. 

Published in Journal of Physics: Conference Series 30 (2006) 30-40, © copyright 
(2006) IOP Publishing Ltd. 



PACS numbers: 02.90.+p, 03.65.Fd 

AMS classification scheme numbers: 16W22, 14L24 

1. Introduction 

In many physical applications one has to do with functions that are invariant under 
transformations of a compact symmetry group. These functions may be considered 
as functions defined on the orbit space of the group and some arguments, like the 
study of symmetry breaking or the study of phase transitions, are better understood 
if studied in the orbit space. To study invariant functions in the orbit space one has 
to know first the equations and inequalities that define the orbit space and its strata. 
A standard way to find these equations and inequalities is by using an integrity basis 
and the P-matrix. This method is reviewed in sections |2] and 0] The orbit spaces 
have not been widely employed in the study of invariant functions yet because the 
integrity basis is often unknown (the determination of an integrity basis of a general 
group is still an open problem of constructive invariant theory) . Section |2| reviews 
some selected arguments of invariant theory and section (the main section of this 
article) shows how the P-matrix can be used together with the Molien function to 
determine the integrity basis of a compact group. 

2. Orbit Spaces 

In this article G is a compact group acting effectively in a finite dimensional space. 
Then, in all generality G is a group of real orthogonal matrices acting on the real 
vector space R™. 
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The orbit through a point x G R n is the set: f2(x) = {g-x, Vg € G}. Orbits do not 
intersect and define a partition of R™. The isotropy subgroup (or stabilizer) G x of the 
point x G R" is the subgroup of G that leaves x unchanged: G x — {g E G \ g ■ x — x}. 
Isotropy subgroups of points in a same orbit are conjugated: G g . x = g ■ G x ■ g~ x . The 
orbit type of the orbit Sl(x) is the conjugacy class [G x ] of its isotropy subgroups. The 
orbit types define a partition of orbits (and of R n ), each equivalence class Ergi is 
called a stratum (of R n ): E[#] = {x E R" | G x € [H]}. The orbits and the strata 
can be partially ordered according to their orbit types. The orbit type [H] is said to 
be smaller than the orbit type [A] : [H] < [A] , if W C K' for some W E [H] and 
A' E [A] . Then [A"] is greater than [H] . The greatest orbit type is [G] and its stratum 
Ergn contains all fixed points of the G-action. Ergn coincides with the origin of R™ if 
the G-action is effective. Due to the compactness of G, the number of different orbit 
types is finite and there is a unique smallest orbit type [G p ] , called the principal orbit 
type. The stratum E[g p ] is called the principal stratum, all other strata are called 
singular. 

The orbit space is the quotient space R n /G. The natural projection tt : R™ — > 
R™/G maps orbits of R" into single points of R™/G. Projections of strata of R" 
define strata of R"/G. The principal stratum of R™/G is always open connected and 
dense in R"/G. If [A"] > [H], then 7r(£rK-i) lies in the boundary of 7r(E[ ff j) and the 
boundary of the principal stratum contains all singular strata. 

A G-invariant function / : R™ — > R is such that f(gx) — f(x), Vg E G , x E R™. 
/ is then constant on the orbits and it is in fact a function defined on the orbit space. 

Some symbols. R[R"] and R[R™] G indicate the rings of real polynomial functions 
and of real G-invariant polynomial functions defined in R", respectively; R[R"]d and 
R[R"]f indicate the rings of homogeneous polynomial functions of degree d in R[R n ] 
and in R[R"] G , respectively; R[pi, . . . ,p q ] indicates the ring of polynomials with real 
coefficients in the q indeterminates pi, ■ ■ ■ ,p q - 

By Hilbert's theorem, R[R™] G has a finite number q of generators pQ: 

VpeR[Rf,3p€R[ Pl ,...,p,] | p(x)=p( Pl (x),..., Pq (x)), VxeR". 

The q generators for R[R™] G are called basic invariants or basic polynomials and they 
form an integrity basis for G (and for R[R™] G ). The integrity basis pi, . . . ,p q is said to 
be minimal (abbreviated in MIB), if no p ai a = 1, . . . , q, is a polynomial in the other 
elements of the basis. The basic invariants p%, . . . ,p q may be supposed homogeneous 
and ordered according to their degrees d a , for example d a < d a +i- One may choose, in 
all generality, pi(x) = Y)j—y xf. The choice of a MIB is not unique, but the numbers 
q and d a ,a = 1, . . . , q are uniquely determined by G. 

A MIB can be used to represent the orbits as points of R 9 . In fact, given an orbit 
ft, the vector map p(x) = (pi{x), P 2{x), . . . ,p q {x)) is constant on fl, and, given two 
orbits fli ^ 0,2, 3 p a G MIB | p a (xi) ^ p a { x 2), Vxi G fti, X2 G ^2- p(x) then defines 
a point p = p(x) G R 9 that can be considered the image in R 9 of fl because no other 
orbit is represented in R 9 by the same point. 

The vector map: p : R™ — > R 9 : x — > ( P i(x), P 2(x), . . . ,p q (x)), is called the orbit map 
and maps R™ onto the subset S = p(R n ) C R 9 . p induces a one to one correspondence 
between R"/G and 5 so that S can be concretely identified with the orbit space of 
the G-action. All the strata of S are images of the strata of R n through the orbit map 
and if [K] > [H], p(T,[ K ]) lies in the boundary of p(E[jjj). The interior of S coincides 
with the image of the principal stratum and the image of all singular strata lie in the 
bordering surface of S. p(E[ G j) always coincides with the origin of R 9 . 
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Generally, the q basic invariants are algebraic dependent, that is, some polynomial 
/ G R[pi, • • ■ ,p q ] exists such that f(pi(x), . . . ,p q (x)) = 0,Vx S R". In this case the 
polynomial f(pi, ■ ■ ■ ,p q ) is called a syzygy (of the first kind). Given a syzygy /, the 
equation f(p±, . . . ,p q ) = defines a surface in R 9 and S must be contained in that 
surface. So, S C Z, where Z is the intersection of all syzygy surfaces. 
All real G-invariant C°°-functions (not only polynomial functions) can be expressed 
as real G°°-functions of the q basic invariants of a MIB and hence define G°° 
functions on R 9 : f(x) = f(pi(x), . . . ,p q (x)) -> f(p) = /(pi, . . . ,p q ). The functions 
f(p) are defined also in points p S but only the restrictions f(p) \ pe s have the same 
range as f(x), x £ R". All G-invariant C°° functions can then be studied in the orbit 
space S but one needs to know all equations and inequalities defining S and its strata. 

Details and proofs of the statements here recalled may be found in (21 BJ 03 El an d 
references therein. 

3. Invariant Theory 

Given a MIB p±, . . . ,p q , there are in general s < q basic polynomials that are algebraic 
independent: Phx, ■ ■ ■ ,Ph,- H s < q necessarily there are syzygies (of the first kind). 
One can determine a finite minimal set of r\ linearly independent homogeneous 
polynomials, . . . , frl' € R[pi, ... ,p q ], of degrees m\\ < ... < m\ Tx with respect 
to the natural grading assigned to the p a (deg(p a ) — deg(p a (x))), that generate all the 
syzygies of the first kind as a module over R[pi, . . . ,p q ]. One can continue to define 
syzygies. Linear combinations of the r± variables f frl\ with coefficients in 
R[p/j 17 . . . ,Ph 3 ], that vanish identically by substituting the with their expressions 
in terms of the p a , are called syzygies of the second kind. Let fi^j-.-jfrf 1 be a 
linearly independent basis for the syzygies of second kind, made up by homogeneous 

polynomials (with respect to the natural grading) of degrees 777,21 < • • • < "i2r 2 • Linear 

(2) (2) 

combinations of the r2 variables /{ , . . . , fr 2 , with coefficients in Rfp/ij , . . . ,Ph s ], that 
vanish identically by substituting the fj^ with their expressions in terms of the f^\ 
are called syzygies of the third kind, and so on. There exist a number m such that one 
has syzygies of the m-th kind f[ , . ■ ■ , /r,™\ but no syzygies of the (777 + l)-th kind. 
Changing the MIB, one changes the syzygies too but does not change the numbers q, 
s, to, d a , Va = 1, . . . , q, rj, to^, V i = 1, . . . ,m, j = 1, . . . ,rj, that are characteristic 



By a theorem of Hochster and Roberts, in R[R n ] G one can choose a maximal 
subset of s algebraically independent polynomials: ph 1 , ■ ■ ■ ,Ph 3 , of degrees d[, . . . ,d' s , 
called primary invariants (or homogeneous system of parameters) , and a number t of 
other homogeneous polynomials: /1, . . . , ft, of degrees toi, . . . ,m t , called secondary 
invariants, such that all p G R[R"] G can be written in a unique way in following 
form: 



where q^ € R[p?!i, ■ ■ ■ ,Ph s ], an d, by definition, /o = 1 [3|Sj- In other words, all 
p £ R[R"] G can be written in a unique way in terms of pt l: ■ ■ ■ ,Ph s , fo, ■ ■ ■ , ft, using 
the fk at most linearly. The most convenient way to choose the ph t and the fk is to take 
Ph t , ■ ■ ■ ,Ph s and fx, . . . , f q -s, equal to the q basic invariants in the MIB, withp;^ = p\, 



of R[R"] 




(1) 



fc=0 
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and /g-s+i, ■■•,/*, equal to some of the products fkfi, with 1 < k < I < q—s, in such a 
way that all products fkfi may be expressed in the form (using in case the syzygies). 
This choice implies that the degrees d[ of the ph t are equal to d a , a G {hi, . . . , h s }, 
and the degrees nik of the fk are equal to d a , a = 1, . . . ,q, a ^ {hi, ■ . . , h s }, or to their 
sums. Then, the degrees mu of the syzygies of the first kind / ; are mu = rrii + mj, 

(2) 

the degrees m-u of the of the syzygies of the second kind /, are m,2i = rrii + rrij +mfc, 
and so on. 

The main problem of constructive invariant theory is to determine a MIB for a 
given group G, and after that all the syzygies. It is then of great interest to find some 
algorithms that suggest what the numbers q, s, d a , m, r^, my, t, m^, are. 

A key role is played by the Molien function M(rf) that may be calculated without 
the knowledge of the MIB in the following way: 

\G\ det(e-ry 5 ) J G det(e - ijg) 

respectively for finite [7] or compact continuous |S] groups, where rj < 1 is an abstract 
variable, e is the unit matrix of G, \G\ is the order of G, dg is the normalized Haar 
measure and the integration is over the compact group variety. The integral may be 
reduced to a sum of integrals over unit circles in the complex plane and calculated 
using the residues ^3 m . An alternative algebraic method to compute the Molien 
function for compact connected groups is presented in |12j . 

When one calculates the Molien function from the formulas written above, one 
ends up with M(rj) expressed by a rational function. 

The expansion of M (rf) in power series has the following interpretation: 

oo 

M 1 (r 7 )=^dim(R[R»]f)r 7 d , 

d=0 

where dim(R[R"] d ? ) is the dimension of R[R"]^ as a linear space. Let's call Mi (rj) 
the first form of the Molien function. Mi(rj) gives a first information about the 
number and degrees of the basic invariants of small degrees. One has dim(R[R"]^) = 
n d°^ ~ n ^ n< d ^ ~ • • • ' wnere n d°^ ls t ne number of linearly independent polynomials 
of degree d, is the number of linearly independent syzygies of the first kind of 
degree d, and so on. For sufficiently small d, one may find out the decomposition of 
dim(R[R™] d ; ) in terms of , • ■ ■ , but in general one is not able to deduce the 

numbers n d l \ 

Two rational forms of the Molien function are particularly interesting. 
The rational function 

Ili=i( 1 -»rO 

here called the second form of the Molien function, has the denominator that is the 
product of s factors of the kind (1 — r\ '), one for each primary invariant ph n and the 
numerator that is a reciprocal polynomial that gives the number t of the secondary 
invariants fk, and their degrees rrik- 
The rational function 

M3(r?) = m =1 d-^) ' 
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here called the third form of the Molien function, has the denominator that is the 
product of q factors of the kind (1 — r] da ), one for each p a in a MIB, and the numerator 
that is a polynomial that gives a partial information about the number of the 

(i) 

syzygies of i-th kind /j , and their degrees my. 

When one calculates M(rj) one ends up with a rational function that is generally 
different from both M<2(rj) and M^{rj). To get M<i{r}) and M^{rj) from M(rf) one has to 
multiply numerator and denominator of M{ti) by convenient factors. Obviously, one 
has M{n) = M x {rj) = M 2 {rf) = M 3 (r?). 

If G is finite then s = n. If G is a compact Lie group then s = n — dim G + dim G v . 

The expressions for Mi(r\) and M^irf) coincide only if q = s (and this implies 
t = m = 0). Many groups for which q = s are classified in 1141 ITS1 . 

In the general case s < q, M2{rf) and M^rf) are unknown as long as one 
does not know the MIB completely. A possible expression for M2(rj) can easily be 
found by requiring only positive coefficients in the numerator. At this point other 
possible expressions for -/^(ry) are obtained by multiplying the numerator and the 
denominator of M2{rj) by factors like (1 + r] d '), where d[ appears in the exponents of 
the denominator, increasing in this way the number of secondary invariants required 
by the numerator. Generally, when one starts to calculate the MIB, one knows the 
number s of algebraically independent polynomials, so one is able to write down the 
correct form for M2{rf). For a given expression of M2(r]), one has to guess the correct 
expression of M 3 (i]). One notes that: 

(i) One obtains Mz{rf) by multiplying numerator and denominator of M2(r)) by a 
certain (unknown) number (q— s) of factors like (1 — rj da ), one for each basic non- 
primary invariant fk, k = 1, . . . , q — s (and this implies the equality m = q — s). 

(ii) The coefficients in the numerator of M 3 (r]) corresponding to the syzygies of the 
first kind are negative. 

(iii) An invariant fk, whose degree vrik appears in the numerator of M2{rj), must 
coincide or with a basic non-primary invariant, or with a product of them. This 
implies that: (i) for each fk, whose degree cannot be written as a sum of 
degrees of basic polynomials, there is a factor (1 — rf ak ) in the denominator of 
M 3 (r]); (ii) the degrees of the syzygies of the first kind are > 2m\. 

If one does not know the numbers q and d a , there are many "possibly right" 
expressions for M^(rf) that satisfy all the points written above, but only one is the 
correct one, corresponding to the right numbers q and d a . Moreover, given the right 
expression for M.${r\)^ the numbers m, r^, my, are not immediately readable from the 
numerator because of simplifications between similar terms. In any case, if there is 
only one kind of syzygies of degree d (for example if d is small) , the coefficient of n d 
in the numerator allows to know the number of independent syzygies of degree d. 

Summarizing, the Molien function is a first powerful device to study R[R"] G 
because it gives some information about the degrees d a and allows, in any case, to 
determine correctly the number and degrees of the basic invariants of small degrees. 

When one knows the degrees of some invariants, one may find them explicitly 
using known properties of the group action or by means of the classical procedure of 
averaging in the group. Given an arbitrary function f(x), its average in the group is 
the invariant function F(x), where: 

F ( x ) = T7^Y, f(3x), or F(x) - / f(gx) dg, 
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respectively for finite |17) or continuous compact |18j groups. Obviously, for the 
linearity of G, if f(x) £ R[R n ]ci, then F(x) £ R[R™]^, so to get an invariant of 
degree d one may take the average in the group of a monomial of degree d. If d and 
n are small, one has only a limited number of possible choices of monomials of degree 
d, but, if d or n are great, this averaging procedure becomes difficult to carry out, 
because of the many possible monomials of degree d one may choose to average and 
because most of the invariants that one finds out in this way are not independent from 
those that one already knows. 

For a detailed presentation of invariant theory and for the proofs that are here 
omitted, the interested reader may read P EE El EH HOI EH an d the references 
therein. 



4. P-matrices and Orbit Spaces 

Given a stratum £ C R™, in a point i£E, the number of linear independent gradients 
of the basic invariants is equal to the dimension of the stratum p(S) C S [221 - O ne 
may calculate the q x q Grammian matrix P(x) with elements P a b{x) that are scalar 
products of the gradients of the basic invariants |23j : 

i— 1 1 1 

so, if x e £, rank(P(ir)) = dim(p(S)). 

From the covariance of the gradients of G-invariant functions (V/(g • x) 
ij ■ V/(x)) and the orthogonality of G, the matrix elements P a f,(x) are G-invariant 
homogeneous polynomial functions of degree d a + db — 2, so they may be expressed as 
polynomials of the basic invariants: 

Pab(x) = P a b(pi(x), ...,p g (x)) = P a b(p) Vie S R" and p =p(x). 

One can define then a matrix P(p) in R 9 , called a P -matrix, having P a b(p) for elements. 
At the point p = p(x) £ S, image in R 9 of the point x £ R" through the orbit map, 
the matrix P{p) is the same as the matrix P(x), P(p) is however defined in all R 9 , 
also outside S, but only in S it reproduces P(x), Vx £ R™. 

P(p) is a real, symmetric qxq matrix and its matrix elements P a b{p) are homogeneous 
polynomial functions of degree d a + db — 2 with respect to the natural grading, 
rank(P(p)) is equal to the the dimension of the stratum containing p, and S is the 
only region of Z where P(p) is positive semidefinite |rjl I24| . 

If s < q basic invariants are algebraic independent, one has rank(P(p)) = s, Vp in the 
principal stratum, and S has dimension s. 

The matrix P(p) completely determines S and its stratification. Defining Sk the 
union of all fc-dimensional strata of S, (As = 1, ... , s), one has: 

S = {p£Z\ P{p) > 0}, S k = { P £Z\ P(p) > 0, rank(P(p)) = k} 

To find out the defining equations and inequalities of fc-dimensional strata, one may 
impose that all the principal minors of P(p) of order greater than k are zero and that 
at least one of those of order k, and all those of order smaller than k are positive. 
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The P-matrix contains all information necessary to determine the geometric 
structure of the orbit space and its strata, so, if one wishes to classify geometrically 
the orbit spaces, it is sufficient to classify the P-matrices. This, in turn, gives a non- 
standard group classification. One discovers so that different groups, no matter if 
finite or continuous, share the same orbit space structure, despite the fact that the 
strata correspond to different orbit types (some examples are in Table X of US]). They 
exist different MlB's (of different linear groups), with the same number of invariants 
and the same degrees, that determine: (i) the same P-matrix (some examples are in 
Table X of US]); (ii) different P-matrices (an example is given by Entries 6 and 25 of 
Table V of [25]). 

Some examples using the P-matrices to study orbit space stratifications, minima 
of invariant polynomials and phase transitions may be found in |2(il 1251 12 7j . 

5. P-matrices and Invariant Theory 

The P-matrix method to determine a M1B is summarized in the following items (i)-(iv) 
and applied in three examples. 

Let pi, . . . ,p q be a set of q independent G- invariant homogeneous polynomials. 

(i) With pi,...,p q one tries to construct the P-matrix. By Euler's theorem on 
homogeneous functions and the standard choice forpi, one has P\ a {p) = Pii(p) = 
2d a Pa, Va = 1, . . . , q. Taking into account the symmetry of P(p), one has then 
to calculate only the matrix elements P a b(p), with 2 < a < b < q. 

(ii) If it is not possible to express P a b{x) as a polynomial in pi{x), . . . ,p q {x), then 
Pab(x) is an invariant of degree d a + db — 2 that is independent from the known 
ones. One has then to add a new basic invariant to the set pi, ■ ■ ■ ,p q . One may 
choose, for example, an irreducible invariant factor of P a b{x) independent from 
Pi,--- ,Pq- 

(iii) When one succeeds to express all the matrix elements P a b(x) in term of the basic 
invariants, one has found a "closed" P-matrix P{p) and a set p\, . . . ,p q that is a 
MIB (or part of it) for a group G that either coincides with G or contains G as 
a subgroup (see Example 3). 

(iv) To check if p\ , . . . , p q do form a complete MIB for G one has to look to the Molien 
function of G. If Ivhivi) and M^rj) agree with the MIB p\,...,p g , then this is 
the MIB for G one was searching in. If instead M3 (77) requires a different set of 
degrees for the basic polynomials, then one has to find out (by averaging in the 
group or in some other way) at least one of the missing basic polynomials of the 
lowest degree. If po is this new invariant, first of all one has to check if the higher 
degree invariants in the set pi, . . . ,p q are all independent from p , otherwise one 
has to discard the non independent ones. Then, return to item (0) above with p 
in the set of the q invariant polynomials. 

Some technical suggestions. With the P-matrix method, with only invariants of degree 
2 one is not able to determine higher degrees invariants and with only even degree 
invariants one is not able to determine odd degree invariants. So it is convenient to 
start to build the P-matrix having determined in some way at least all invariants of 
degree 2, all invariants of the smallest odd degree and all the invariants of the smallest 
even degree greater than 2. One then starts to determine the P-matrix as indicated 
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in items <{i|)- (|rv|l above, proceeding step after step to calculate the matrix elements of 
increasing degrees. During this calculation, every time one adds a basic invariant to 
the set of basic invariants p%, .. . ,p q , according to item (|n| above, it is convenient to 
look if this MIB is consistent with the guessed expression of M2(r/) and Mz(rj). In 
case, one has to modify the expressions of M 2 (?7) and/or M 3 (r)) to include the new 
basic invariant. One has so a partial check of the correctness either of -^2(77) and 
AI3 (77) , either of the MIB. Every time one starts to determine the matrix elements 
Pab(p) of a given degree d, it is convenient to control if one knows explicitly all the 
basic invariants of degree lower than d that are predicted by M^r/). If one misses one 
of these invariants it is strongly recommended to find out this basic invariant first, 
otherwise, one does longer and harder calculations. 

If a subset pi ± , ■ ■ ■ ,pi k of the MIB for G is a complete MIB for a group G C 
0(n, R), containing G as a subgroup, then the principal minor of the P-matrix of G 
formed by the rows and columns l\, . . . ,lk is expressed only in terms of pi t , . . . ,pi h . 



Example 1. Let G — Oh, the symmetry group of the cube. The Molien function of G 
has the following expression: 



M 2 {rj) = M 3 (rj) 



1 



(1 - 77 2 )(1 - 77 4 )(1 - 7? 6 ) 



and predicts a MIB with 3 algebraic independent polynomials of degree 2, 4 and 6 
that have been calculated by many authors (see for ex. Table 4 of As an elemen- 
tary example of the P-matrix method to determine a MIB, let us start with the basic 
invariants pi — x 2 + y 2 + z 2 and P2 = x 4 + y A + z A . The only non trivial element of the 
P-matrix formed with p\ and P2 is P 2 2- One finds easily Vp2 = (4a; 3 , 4j/ 3 , 4z 3 ) and 
P22 = Vp2- Vp 2 = 16(x 6 + y 6 + z 6 ). As P22 cannot be written in the form cip\+C2PiP2, 
the basis pi and P2 is not complete and one may take p^ — x 6 + y 6 + z 6 for a third 
element of the MIB. All elements of the P-matrix constructed with p\, P2 and p 3 can 
be written only in terms of p\, p 2 and P3, so the MIB is complete, it agrees with the 
Molien function and in fact is a possible MIB for G. This trivial Example shows that 
by constructing the P-matrix with only part of the MIB known, one succeeds to find 
out all the missing basic invariants. 



Example 2. Let G be the group of order 192 generated by the matrices 



,9i 



This group is the realification of the complex unitary group considered in and 
one is referred to f° r details about its structure. 

The Molien function of G has the following first and second form: 





/ 1 


1 





^ 




( 1 








\ 




1 


-1 



















-1 


2 








1 


1 


92 = 








1 







V 





1 


- 1 J 




V 


1 





/ 



Mi (?y) = 1 + if + rf + rf + 4t? 8 + 4r/ 10 + 7r] 12 + 7if 



12n 



10 



-1777 18 + 22?f u + 22r] 22 + 36i] 2 



M 2 (v) 



1 + 2^ + 2r] LZ + 2rj Lb + 5rj 



5r) 2 



- 43t7 2 ' 
2ry 26 



2t7 3U + 2^ 



(1 - 77 2 )(1 - 77 8 )(1 - r/ 12 )(l - 77 24 ) 
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From M\(n) one gets the following kind of information. The term Aif suggests that 
there are 4 lin. indep. invariants of degree 8: p\, P2, P3, Pa- The term 7r] 12 suggests 
that there are 7 lin. indep. invariants of degree 12: p\, p\p2, P1P3, PiPa, P5, Pe, 
p-j. The term 12n 16 suggests that there are 12 lin. indep. invariants of degree 
16: 13 lin. indep. invariants are the terms in the expansion of the polynomial 
(P2 + P3 + Pi + p\) 2 + {P5 + P6 + P7)pI and it must exist one syzygy of the first 
kind of degree 16, so that 12 = 13 — 1. The term 17r/ 18 suggests that there 

are 17 lin. indep. invariants of degree 18: 12 independent invariants are found by 
multiplying by p\ the 12 independent invariants of degree 16, and there must exist 
5 new basic invariants of degree 18: p%, p$, pio, Pn, P12, so that 17 = 12 + 5. One 
may continue this analysis and find out that there are 4 syzygies of the first kind of 
degree 20: ■ ■ ■ , and no basic invariants and no syzygies of degree 22. The 
analysis becomes very uncertain for degree 24 because one might have either basic 
invariants, either syzygies of the first and of the second kind, and one only knows that 
their numbers n 2 °4 , n^} and n 2 2 } must satisfy the condition n 2 °^ — n^} + fi^l = 36. 

The information obtained so far from M\ (rf) and M 2 (rj) suggest to obtain the third 
form of the Molien function by multiplying numerator and denominator of M 2 (n) by 
(1 — jf ) 2 (l — j] 12 ) 2 (I — rj 18 ) 5 . Note that the exponents 2, 2, 5 are the coefficients of 77 s , 
n 12 , n 18 in the numerator of the second form of the Molien function. Doing this, one 
finds a third form of the Molien function that shows at the numerator the existence of 
the syzygies f^, . . . , but suggests, wrongly, that there is only one basic invariant 
of degree 24. At this point it is convenient to stop the (uncertain) analysis of the 
Molien function and start to build the P-matrix. 

The Molien function suggests so far a MIB with 1 invariant of degree 2: pi, 3 
invariants of degree 8: p2,P3,P4, 3 invariants of degree 12: P5,Pe,P7, 5 invariants of 
degree 18: P8,P9,Pio,Pu,Pi2, and one invariant of degree 24: pi 3 . In this example 
the explicit calculation of all the invariants by the method of averaging in the group 
is not trivial, because of the high degree of the invariants, so, it is convenient to use 
the P-matrix method. 

First of all, one defines: pi(x) = x 2 +x 2 +x 2 +x 2 and finds out the three invariants 
of degree 8, for example by averaging in the group the monomials x\x 2 x%x\, x\x\ and 
x\x\. One then starts to build up the P-matrix. One calculates first the matrix 
elements of degree 14: P2,2(x), ■ ■ ■ , Pa^(x) and one finds out the 3 basic invariants of 
degree 12. One then calculates the matrix elements of degree 18: P2,s(x), . . . , Pt^x), 
and one finds out the 5 basic invariants of degree 18. One then calculates the matrix 
elements of degree 22: P^^{x) 1 . . . , Pj^(x), with no need of new basic invariants. One 
then calculates the matrix elements of degree 24: P2,s(x), . . . , Pt,i2(a;), and one finds 
out two basic invariants of degree 24: ^13(2;) and pu{x), not only one, as expected. 
Then the third form of the Molien function has to be corrected by multiplying its 
numerator and denominator by (1 — r/ 2 ) and obtaining: 



1 — r\ 



,16 



4?? 20 _ ^24 _ 8?? 26 + 4r? 28 _ gr/ 30 + g?? 34 _ 



19?f 6 + . . . 



M 3 (v) = 



(1 - r7 2 )(l - ?7 8 ) 3 (1 - 77 12 ) 3 (1 - 77 18 ) 5 (1 - t? 24 ) 2 



Calculating all the remaining elements of the P-matrix one does not need any 
new basic invariant. As pi, . . . ,pu agrees with the third form of the Molien function, 
they form a complete MIB for G. It is important to note that the invariant pu was not 
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predicted from the Molien function, but it is required to express some of the elements 
of the matrix P(x) in terms of the MIB. 

This example shows that the Molien function is not sufficient to specify the 
number and degrees of all the basic invariants in a MIB, but its combined use with the 
construction procedure of the P-matrix allows one to do that. Moreover, the P-matrix 
method is constructive, in the sense that it allows to determine concretely the MIB. 

In this case no proper subset of the MIB pi, . . . ,pi4, is a complete MIB for a 
group G containing G as a subgroup, because no principal minor of the P-matrix can 
be written only in terms of the corresponding subset of basic invariants. 

Having found the MIB, all the syzygies of any kind can be determined. 

Example 3. Let G be the cyclic group of order 3, generated by the matrix 




The three forms of the Molien function are easily found. The second form is 



M 2 (n) 



(i_^)(i-^) 



and suggests a MIB of 3 elements of degree 2, 3, 3. Using the standard form of p\ and 
averaging the monomials xf and x 2 , one finds out the following MIB: 

p 1= x\ + x\, p 2 = x\- 3x x x\, p z =x\- 3x 2 xj . 

If one uses the P-matrix method in this trivial example starting with the basic 
invariants p\ and p 2 , one does not need any other invariant to write down the P- 
matrix and would not find p 3 . The group that leaves p\ and p 2 invariant however 
is G, that has 3 invariants, according to its Molien function. To use the P-matrix 
method it is then necessary to start with all invariants of the smaller degree greater 
than 2. 

Let now H be the cyclic group of order 3, generated by the 4x4 matrix 

9 e 
e 9 

where e is the unit matrix of order 2 and g is the matrix written above. The three 
forms of the Molien function are easily found. The second form is 

1 + 2 V 2 + 6r? 3 + 2 V 4 + ry 6 

2 ™ (1 - ^2)2(1 _ ^3)2 > 

and suggests a MIB of 12 elements: p\,...,pi of degree 2, and p$, . . . ,p\ 2 of degree 3. 
All the basic invariants can be easily found by averaging in the group. For example, 
one may use the standard form of p\ and the averages of the monomials x\, X1X3, 
X1X4, x\, x 2 , £3, £4, X1X2X3, X1X2X4, X1X3X4, X2X3X4 for the invariants p 2 , . . . ,pi 2 . 
One verifies easily that none of them can be expressed as a polynomial in the other 
ones, so they form a MIB for H . 

In this case if one uses the P-matrix method starting with pi , . . . , p 4 and one 
arbitrary invariant of degree 3 one would also find out all the basic polynomials of 
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a MIB (generally with more complicated expressions than those obtained from the 
averages) . 

If one uses the P-matrix method starting for example with the basic invariants 
Pi = x\ +X2+x\ + x\, P2= x\+ x\ and p 5 = x\ — 3xix|, one does not need any other 
invariant to write down the P-matrix. One knows from the Molien function that the 
polynomials pi,P2,P5 do not form a MIB for H and one understands that one has to 
use more invariants at the beginning to find out the MIB with the P-matrix method. 
The invariance group of pi,P2,P5 is the group H, containing H as a subgroup, 
generated by an arbitrary rotation in the plane £3,2:4 and by the matrix g acting 
in the plane x\, x-i- We know from above, or from the Molien function of H: 

M2iV) = (1-^(1-, 3) ' 

that a MIB for H has 2 invariants of degree 2 and 2 of degree 3, so p\,pi,Ph is not a 
complete MIB even for H. 

This example suggests that it is convenient to start using the P-matrix method 
after having determined all invariants of degree 2 and all invariants of the smaller 
degree greater than 2 predicted by the Molien function. Moreover, it shows that the 
Molien function plays a fundamental role to decide if a given set of invariants do form 
a complete MIB for some group or not. 
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